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Abstract. - We derive the transmission coefficient, T(oj), for grazing incidence of crystals with 
spatial dispersion accounting for the excitation of multiple modes with different wave vectors 
k for a given frequency u. The generalization of the Fresnel formulas contains the refraction 
indices of these modes as determined by the dielectric function e(w, k). Near frequencies w e , 
where the group velocity vanishes, T(uu) depends also on an additional parameter determined 
by the crystal microstructure. The transmission T is significantly suppressed, if one of the 
excited modes is decaying into the crystal. We derive these features microscopically for the 
Josephson plasma resonance in layered superconductors. 



Usually propagation of light in crystals is sensitive to average properties described by the 
dielectric function, but not to the specific atomic structure of the crystal, because the wave 
length of light is much larger than the interatomic distance. However, in crystals with spatial 
dispersion the atomic structure may affect optical properties at special extremal frequencies, 
w e , where the group velocity, v g = <9w(k)/<9k, of an eigenmode with dispersion w(k) vanishes. 
Near uj e the effective wave length, X g = v g /uj, becomes comparable with the interatomic dis- 
tance, and then optical properties such as the reflectivity may feel the crystal microstructure. 

Pekar |l| and Ginzburg 0, already realized that several eigenmodes with different wave 
vectors, k, at given u> may be excited by the incident light inside a crystal with spatial disper- 
sion. Then the Maxwell boundary conditions (continuity of the components of the electric, E, 
and of the magnetic field, H, parallel to the surface), are not sufficient to find the relative am- 
plitudes of these modes and hence the reflection coefficient. To overcome this problem within 
the macroscopic approach Pekar and Ginzburg introduced so called additional boundary con- 
ditions (ABC), which are supposed to be related somehow to the crystal microstructure near 
the surface. In this phenomenological approach the choice of these ABC is only heuristic 
and may be controversial, see ref. || and Comments to this paper. It is only a microscopic 
model which can determine the solutions inside the crystal unambiguously and demonstrate 
explicitly how the reflectivity depends on the crystal microstructure. 

In this paper we present for the first time general macroscopic expressions for the transmis- 
sion coefficient T(u) in uniaxial crystals in the vicinity of co e to show how the Fresnel formula 
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is modified by the presence of multiple solutions. Then we find a full microscopic description 
of optical properties for the Josephson plasma resonance (JPR) in layered superconductors, 
which supports the phenomenological picture and derives the ABC. We show explicitly that 
the maximum T(uj) reached at tu ss u e depends on the interlayer distance of the superconduc- 
tor. Weak dissipation as well as a perfect crystal structure are necessary to observe this effect 
of spatial dispersion, as otherwise the influence of the microstructure is smeared out and the 
conventional Fresnel results are restored. For details of the calculation see Q|. 

We also point out that the stopping of light, i.e. v 9 = 0, which in the case of gaseous 
media attracted considerable interest for the coherent optical information storage recently || , 
here naturally appears in a solid at a finite wave vector due to spatial dispersion. 

The JPR is an interlayer charge oscillation due to tunneling Cooper pairs and quasi- 
particles in layered superconductors. It is strongly underdamped because the quasiparticles 
responsible for dissipation are frozen out at low temperatures. The JPR is an appropriate 
phenomenon to study the effect of spatial dispersion on optical properties both theoretically 
and experimentally. Firstly the quite simple Lawrence-Doniach model [|| Q is sufficient to 
provide a complete microscopic description. Secondly, weak dissipation may be achieved at 
low temperatures if the JPR frequency is well below the superconducting gap. Thirdly, there 
is experimental evidence that the spatial dispersion of the Josephson plasmon is significant. 
Recently in the layered superconductor SmLai_ K Sr :E Cu04_a two peaks at sa 7 and rs 12 cm -1 
were observed in optical properties H . These peaks can be naturally understood as the JPR of 
alternating intrinsic junctions with SmO or La^^O^ in the barriers between the CuCvlayers, 
and the peak intensities indicate a strong inter-junction coupling P,[^,[s|-[Tl|], i.e. a significant 
spatial dispersion of the plasma modes. Tachiki et al. mentioned the excitations of multiple 
modes and Marel et al. suggested an effective dielectric function for alternating junctions, 
but in neither case the effects on the reflectivity presented below were discussed. 

To formulate the general problem posed by spatial dispersion on a macroscopic level we 
consider a layered crystal with the dielectric functions e a (w) along the layers and e c (u),k z ) 
along the z-axis. In e c (u>, k z ) we account for a collective mode (here JPR) with the dispersion 
ui c (k z ), i.e. e c [(jj c (k z ), k z ] = 0. The incident polarized light of frequency ui has a magnetic field 
B = (0, B y ,0) and a wave vector ko = (uj sin 9/c, 0,ui cos 6/c), where 9 is the angle relative to 
the z-axis normal to the surface and to the layers. The bulk Maxwell equations 

ck x B y = -uje c (uj, k z )E z , ck z B y = ue a (oj)E x , k x E z - k z E x = -(w/c)B y , (1) 

determine the wave vector k = (ujsm9/c, 0, k z ) of the eigenmodes at the given frequency u>. 
Here k z {ui) is a solution of the equation 

k 2 z = e a H( W 2 /c 2 )[l - sin 2 e/e c (u, k z )}. (2) 

When cu c and e c are /^-independent, eq. @ has one solution for k\ and the Maxwell boundary 
conditions lead to the conventional Fresnel formula T = 1 — |1 — k| 2 /|1 + k\ 2 with k = 
n/[e a {u>) cos 9], where n = ck z (u))/u> is the refraction index. Then T(u>) is peaked near u) c . 

In a crystal with spatial dispersion in e c ((J,k z ), Eq. (|J) has multiple solutions for k z (uj), 
in the simplest case four (real or complex) ones, ±ni(w), ±77,2(0;) Jl|, |), which can interfere 
with each other ||] . For a semi-infinite crystal at z > only two of them are physical. At low 
dissipation the energy flow (Poynting vector S) of proper modes should be directed into the 
crystal in order to preserve causality. Thus their group velocity, v gz (uj) — duj(k z ) / dk z ~ S z , 
should be positive 0. In the case of normal (anomalous) dispersion this requires positive 
(negative) Re(fc z ). When dissipation is taken into account, this rule is equivalent to the 
condition that the eigenmodes decay inside the crystal, Im(fc z ) > 0. Without dissipation one 
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gets v gz = at the frequency io e , separating the normal and anomalous parts of the dispersion 
curve uj(k z ). At lo — lo e , the two physical solutions for k z have the same \k z \, but opposite 
signs, i.e. ni(u> e ) + n2(w e ) = 0. 

To determine the relative amplitudes of the two modes excited by incident light, ABC 
for the macroscopic polarization P were used in In the bulk P z is related linearly to 

E z by the susceptibility Xc(k z ) — [e c (w, k z ) — e c o]/47r, e c o being the high frequency dielectric 
constant. Generally this relation may differ near the surface. This difference is very small for 
the JPR, as it is confined between adjacent superconducting layers both in the bulk and on the 
surface, and the surface intrinsic junction is almost the same as that in the bulk. Therefore 
we restrict ourselves to the dielectric functions e c (z,z') — Q(z)Q(z')e c (z — z') with a sharp 
cutoff at the surface z = 0. 

The ABC for our system, as proposed by Ginzburg, 

P z (z)+£(dP z /dz) = 0, z^O, (3) 

determine the amplitudes 7^2 of the fields 

E z (z) = 71 exp(ifc z:L z) + j 2 exp(ik z2 z), (4) 
Pz(z) = 7iXc(fc z i)exp(i/c z iz) + ~f 2 Xc(k z2 ) ex.p(ik z2 z). (5) 

Here I is a phenomenological parameter related to the crystal microstructure and is to be 
determined in a microscopic model. Usually, Xc(k z ) — Xc(0) oc k z at k z — > 0. In this limit the 
relation for the amplitudes is 

7i(l + + 72(1 + i£n 2 ) = 0, i = uijc. (6) 

Using eqs. (0)-(|)) and @, we derive 

1 \ + e- x n x n % 

K = 

cost) m + n 2 ~ i£,e a (l + e a nin 2 ) 

In the limit \n 2 \ 2> \ni \ 3> -\/|e a | the refraction index with smallest \n p \ determines k and this 
condition defines the frequencies where the one-mode Fresnel description is correct. In the 
two-mode frequency interval T(uj) is peaked when \n\ \ — \n 2 \, and two situations are possible. 



When both rc.1,2 are real without dissipation we have a frequency u> e , where n\ + n 2 = 0. 
At w = w e only the term i^n\n 2 with £ <c 1 remains in the denominator leading to T(uj e ) — 0. 
As lo increases above ui e , the transmissivity increases and reaches its maximum, 

T max = T(w e , max ) = 2/[(l + e 2 a e cos 2 B) 1 ' 2 + 1], (8) 



at lo — w e ,max when n\ + n 2 = (1 + e a ; 1 rjin2)(cos _2 8 + e 2 ^ 2 ) 1 / 2 above u) e . We see that T max 
depends on £ and is generally smaller than the Fresnel result T max = 1 without dissipation. 
Note, that the opposite signs of the refraction indices ni :2 near Lo e due to causality are essential 
for the dependence of T max on the length I. The vanishing of n\ + n 2 at to e in eq. ([7]) and its 
consequence on T max have not been realized previously [l2| . 

We anticipate an even stronger suppression of T max due to spatial dispersion at frequencies 
LOi near lo c , where n\ is real, while n 2 = in\ without dissipation. This occurs in supercon- 
ductors when the dispersion of the collective mode is anomalous. T(oj) is peaked at cj,, 
but T(oji) <C T(Lo e , mayi ) (eq.(7)), as now m + n 2 ^ in the denominator of eq. (0). The 
microstructure ~ ^ is irrelevant at Wj. This observation is confirmed below for the JPR. 



(7) 
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The frequencies ui e and u>i, the frequency interval of the two- mode regime near u> e and w%, 
as well as 71^2 are determined by the dielectric function, but T max near Lo e and oji also depend 
on the ABC. These features and the ABC will now be derived microscopically for the JPR. 

This requires the solution of the Maxwell equations and the equations for the phase of 
the superconducting order parameter, exactly accounting for the atomic (layered) structure 
of the crystal along the z-axis, supercurrents inside the 2D layers at z = ms and for interlayer 
Josephson and quasiparticle currents between neighboring layers. We obtain 



cd z B v = ie a0 uj 



E x ~^fY E x sS(z - ms) 

m=0 



d z E x — ik x E z = i—B y ,E zm , m+ i 
c 



(m+l)s dz 

E z — . 

s 



ck x By — 



E z 



m=0 



mj m,m+l 



Cjfes 



m,m+l 



Pr, 



(9) 
(10) 

(11) 
(12) 



Here ji m is the chemical potential in the layer to, V m>m+ i is the difference of the electro- 
chemical potentials, uj^ = 87r 2 csJ;/e c o are the bare JPR plasma frequencies determined by 
the Josephson critical current densities J; in junctions of type I = 1,2, uj a o — c/ X a by/e a o 
is the in-plane plasma frequency, e a Q is the high frequency in-plane dielectric constant and 
X a b is the in-plane London penetration length. The function f i s defined as fm,m+i{z) = 1 
at ms < z < (to + 1)s and zero elsewhere. To obtain eq. (hll) for small amplitude oscil- 
lations we expressed the supercurrent density J„ ra+1 = Ji sin (p m ^ n +i sw Jnp m ,m+i via the 
phase difference <p m ,m+i — ^iV m ,m+i/^- Further, ujf = lo 2 {1 — iAira c iui / w^e^) -1 takes into 
account the dissipation due to quasiparticle tunneling currents, J^m+i = cr ciVm,m+i/ es j de- 
termined by the conductivities a c i. We express the difference p m — p m +i via the difference 
of the 2D charge densities, p m , as p m — fx m +i — (47rsa/e c o)(p m — p m +i), where the param- 
eter a = (e c o/4'Kes)(dfj,/dp) characterizes the inter-junction coupling, i.e. the JPR charge 
dispersion ||. This parameter a is estimated as 0.4 for SmLai-zSr^CuO.}-.? |10|. Finally, the 
Poisson equation expresses p m via E z . 

The solution of eqs.(|9j) - (12) between the layers to and to + 1 is 

E z {z) = c m exp(igz) + d m exp(-zgz) + aA m , (13) 

and similar for B y and E x , where g = (w 2 e a o/c 2 a) 1 / 2 and a -1 = 1 — sin 2 9/e c o- The continuity 
of By and E x across the layers leads to a set of finite-difference equations for c m , d m and A m . 

First we consider a crystal with identical junctions, luqi = ujq = c/X cy /e c o and <j c i = crc- 
Omitting the terms of order (gs) 2 « (s/A c ) 2 ~ 10~ 10 we obtain in the bulk 



A m (w - a - 2a) + a{A m +i + A r 
at to > 1 and similar equations for c m 



_i) = (l-2a/3)(c m + d m ), (14) 
d m . Here w — Co 2 /lo 2 an d ft = s 2 l^a\b a 



for A„ 

(~ 10~ 4 in cuprates). Eliminating c m , d m we find linear equations for A m alone, which give 
the dispersion of the eigenmodes. The difference between the equations for A and Ai is 
the microscopic boundary condition which allows us to find A m using bulk equations. The 
dispersion of the eigenmodes (0 < q < 2ir), 

L0 2 (q) 



w(q) = 



= 1 + 2a(l - cosg) 



(a -1)13 
(3 + 1 — cos q ' 



(15) 
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Fig. 1 - Schematic picture of the dispersion w — lu 2 /uiq for two alternating junctions (a c — 0, solid 
line). The lower band is analogous to the case of identical junctions and its dispersion is normal for 
8 — (dashed). Its mixing with a decaying electromagnetic wave (as shown by the dashed, vertical 
line at sin 2 (g/2) = —2/3) results in two propagating modes with real q near the lower band edge w e , 
where the group velocity vanishes. The anomalous dispersion in the upper band gives rise to one 
propagating and one decaying mode and a special point Wi, where i%2 = mi. 

is equivalent to eq. (||) with 

e c (u>,q) = e c0 [l - uj 2 c (q)/u 2 + 47rier c /we c0 ] , e a (oj) = e a0 [l - u 2 a0 /uj 2 ] , (16) 

where uj 2 (q) = ui 2 [l + as 2 k 2 ] and k 2 = 2(1 — cosq)/s 2 . The dispersion of the plasma mode 
due to the charging effect, u c (q), is normal at low q Q. The last term in eq. fll5| ) describes the 
anomalous (at small q) dispersion due to the inductive coupling of the in-plane currents |tJ . 
When alone, it describes the decay of an electromagnetic wave inside the superconductor 
on the length X a i,. As shown in fig. 1 for the lower band, their combination leads to the 
extremal frequency, w e = 1 + u, with u — [8 (a — l^a] 1 / 2 , if a > (a — l)/3/8 and if the 
dissipation is weak, 47r<7 c /woecO <C u. For any angle 9^0 the longitudinal plasma oscillations 
mix effectively with the electromagnetic wave at small q \JuJ1a. From eq. ( [Tq) we find 
n i,2 = (c/^s) 2 [w - 1 ± y/(w - l) 2 - u 2 ]/2a and \nxn 2 \ = A 2 e c0 u/2as 2 w A 2 /(sA q6 ) > 1. 
Hence, the two-mode regime, when |m| ~ \ n 2\, occurs in the frequency interval of width 
w uloo above u> e . Otherwise, at \w — 1| ^> u, the standard Fresnel formulas are correct. 

Note that by modifying the JPR frequency tuo and hence ui e , e.g. by an external magnetic 
field, the stopping of light at lu c can be modified fastly, which might serve as the building 
block of a future magnetooptical device. 

Next we find solutions for A m , c m , d m . We obtain for m > 1 

An = 7i^%i) exp(iqim) + 72^(92) exp(iq 2 m). (17) 

The microscopic boundary condition found as described above is A—± — 0, where A—± is the 
extension of the bulk solution, eq. (p7|), to m = — 1. This corresponds simply to the fact that 
for the surface junction, m = 0, one neighboring junction, m = — 1, is missing. We note that 
A m plays the role of the polarization, P z , because it describes the response of the system to 
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Fig. 2 - Schematic transmission coefficient T(w) with (solid) and without (dashed) spatial dispersion 
(ct c = 0). As a increases, the peak in T(w) is shifted to higher frequency and its amplitude decreases 
depending on the interlayer spacing ~ £ (w = uj 2 /uq, a charge coupling). 



the electrochemical potential, see eq. With A_x = P z (z = — s) we derive microscopically 
eq. (||) with I — —s. Then in eq. we get £e a — sX c ^/e^o/X 2 b which may be of order unity 
in cuprates like Tl-2212 with X c /X a b ~ 100 and the JPR frequency ~ 20 cm -1 . 

Hence, for grazing incidence spatial dispersion (a ^ 0) leads to (i) a shift of the peak 
position in T(w) by u ~ y/]3, and (ii) a decrease of T max depending on the interlayer spacing 
- £, see fig. 2. 

In a crystal with alternating Josephson junctions, 6 — lo^/uj^ < lj we introduce -Ai(n) = 
and A 2 {n) = A^n+i- They describe the polarization inside two different junctions in 
the unit cell n. For a plasmon decoupled from the electromagnetic field (at 8 — 0), the 
charge coupling of junctions leads to the dispersion of the two bands and a gap between them, 
see fig. 1. In the lower (upper) band the dispersion is normal (anomalous). Mixing with 
the decaying electromagnetic mode in the lower band leads to the existence of an extremal 
frequency u> e for a 2 > (a — 1)/3(1 — <5)/4 and 47rcroiA<-'oie c o <C u. Hence, the situation in the 
lower band is similar to that in the case of identical junctions. In contrast to this, the mixing 
of the upper plasma band with anomalous dispersion and the decaying electromagnetic mode 
leads to an eigenmode with anomalous dispersion everywhere. Here rii is real, while rii is 
imaginary, and a critical frequency u>i exists, see fig. 1. 

When a junction of type 1 is at the boundary, the missing neighboring junction is of type 
2 and the microscopic boundary condition is derived as ^(n = —1) = 0. The solution of the 
equations for A\{n) and ^(n), i.e. the generalization of eqs. (14) and ( p7| ) for alternating 



junctions, are the eigenvectors [Ai(q p ), A2(q p )], (p = 1,2), which describe the microstructure 
of the local polarization within the unit cell for each mode p. Here we obtain 

„_ 1 [l + ta ln i n 2][A~2(qi)n 2 - A 2 (q 2 )ni\_ ^ 



cos6 A 2 {qi)n 2 2 - A 2 (q 2 )n 2 1 

where A 2 {q) = A 2 {q) exp(— iq) ps A 2 (Q){1 — iq/2) at small q — 2sk z . This leads in the lower 
band to eq. ^ for k and a similar behavior of T(w) as in the case of identical junctions. In 
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the upper band the atomic scale £ is irrelevant. There the maximum value of T(w) is reached 
at ui max = Wi when ri\ — ini for ctq/ = 0, 



9\ 3 / 2 ^ 



A c (se c o) 1/2 cos( 



(a - l)L 



8a 2 a 



1/4 

(19) 



Here L = Wi(l + S)~ 2 — 8a, Wi — (l + 5)(l + 2a)(l + y/1 ~ p)/2S is the upper edge of the band 
at a = 1, q = and p = 4(5(1 + 4a)/(l + <5) 2 (1 + 2a) 2 . Thus T max depends on a explicitly and 
is smaller by the factor (s/Xab) 1 ^ 2 than T max in the lower band or for the case with identical 
junctions, respectively. This difference vanishes in the Fresnel limit, when 47T(To; ^> ujqiu. 

In conclusion, multiple eigenmodes, propagating or decaying, are excited by incident light 
in crystals with spatial dispersion. At low dissipation this leads to a drop of the maximum 
transmission near resonance frequencies in comparison with the single- mode Fresnel result. 
In the case of two propagating modes near frequencies u> e , where the group velocity vanishes, 
the maximum transmission coefficient, T max , depends on the crystal microstructure, which is 
in remarkable contrast to the general belief that atomic scales cannot affect optical properties 
due to the large wavelength of light. The drop of T max is larger when one excited mode 
is propagating, while the other is decaying. This behavior was demonstrated explicitly for 
the JPR, but it is general for any optically active excitations, e.g. optical phonons with 
anomalous dispersion in insulators. However, the condition of weak dissipation and a perfect 
crystal structure are crucial to observe deviations from the Fresnel regime. 
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